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Abstract 



The time-dependence of correlation functions under the influence of classi- 
rS I cal equations of motion is described by an exact evolution equation. For con- 

^ servative systems thermodynamic equilibrium is a fixed point of these equa- 

tions. We show that this fixed point is not universally stable, since infinitely 
many conserved correlation functions obstruct the approach to equilibrium. 
Equilibrium can therefore be reached at most for suitably averaged quantities 
or for subsystems, similar to quantum statistics. The classical time evolution 
of correlation functions shows many dynamical features of quantum mechan- 
ics. 



^This work was performed in part at ITP, UCSB, Santa Barbara, and supported by the National 
Science Foundation under Grant No. PHY94-07 197. 



Recent years have seen large-scale computer simulations of the statistical be- 
haviour of nonlinear classical field equations. Examples in particle physics and 
cosmology range from the formation of defects during phase transitions in the early 
universe and inflation |^ to the computation of the rate of baryon-number violating 
processes at high temperature T |0 or speculations about the formation of a disor- 
dered chiral condensate in heavy ion collisions [Q. A typical model is a ip'^ theory 
for A^ complex scalar fields, A^ = 1 being relevant for string formation and N = 2 
for the chiral phase transition. Simulations of statistical properties of such systems 
solve the classical nonlinear field equations numerically for given initial conditions. 
Subsequently the average over a suitable ensemble of initial conditions is taken. 
This determines the time-dependent correlation functions of the system. The use of 
classical field equations is often motivated by the conjecture that classical equations 
govern the behaviour of long-distance modes (typical momenta smaller than T) even 
in case of high temperature quantum field theory. We emphasize that quantum ef- 
fects can be (partially) incorporated by using the field equations derived from the 
quantum-effective action or a coarse-grained version of it [Q . 

As a step towards an analytical treatment of the time evolution of statistical 
properties an exact evolution equation for generating functionals of time-dependent 
("equal time"-) correlation functions has recently been proposed ||^. For conser- 
vative classical field equations thermodynamic equilibrium is a fixed point of this 
evolution equation. In this note we investigate if and how this fixed point is ap- 
proached. Intuitively, one expects that suitably averaged "macroscopic" quantities 
should thermalize even in the absence of dissipation: Two types of gas should mix 
due to classical motion with elastic scattering of the molecules.0 Nevertheless, our 
investigations show that the equilibrium fixed point is not universally attractive. We 
find that fluctuations around this fixed point are characterized by infinitely many 
conserved correlations which obstruct the approach to the fixed point. The situation 
seems very similar to quantum statistics where the trace of arbitrary powers of the 
density matrix is conserved for a closed system. Under such cirucumstances only 
suitable "macroscopic" degrees of freedom or particular subsystems can thermal- 

^In case of a finite number of conserved quantities beyond energy the equilibrium should, of 
course, be formulated in accordance with the conservation laws - in analogy to a conserved particle 
number in classical statistics. 



ize whereas on a "microscopic" scale correlation functions fluctuate in time. Our 
treatment of the time evolution of correlation functions in classical statistics indeed 
reveals many analogies to quantum statistics. Our considerations are quite general, 
with applications far beyond the specific (p'^-modeh mentioned in this note. 

Consider some number of complex degrees of freedom Xm whose time evolution 
is given by the equation of motion 

(Summation over repeated indices is always implied.) The matrix H is hermitean 
and we assume that it can be expanded in powers of x- For the purpose of this note 
we will also specialize to the case of a conserved charge (corresponding to phase 
rotations of x) for which H contains equal powers of x and x* 

oo 

J^mn = 2^ '^mnqiri...qiriXqiXri---XqiXri [■^) 

1=0 

Here h^'-^ is symmetric in {m,qi,...,qi) and (n,ri, ...,r;) and obeys h!^mnq^ri...qiri — 
^nLriqi...riqr The equatiou of motion (|I|) can be derived from a quantity U 

oo -| 

^ ~ / . ; I T '''mnqirT...a,r,XmXnXqTXri---Xq,Xri ['j) 

as 

• __-^ - --^ m 

''"'' 'drj ""--'dxm ^^ 

We note that the system is conservative {dU/dt = 0) if the coefficients /i*^'^ are time- 
independent. Eqs. (^, @ are characteristic for non-relativistic field equations. As 
an example we start from the relativistic wave equation for a complex scalar field ip 

dW 

p{x) = {A - M'Mx) - ^{xMx) (5) 

Here VT is a real function of p = f*(f. In the nonrelativistic limit we use 

p{x,t) = e-''''x{x,t) (6) 

and neglect terms ~ x such that 

A 1 dW 
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We identify the infinitely many Xm with x(x) or with corresponding Fourier modes 
x{p) = j cfixe'^^^^'-xix). Obviously H = — ^ + ^^{x) is hermitean. We have 
used a notation where eq. (0) ressembles the Schrodinger equation. It should be 
emphasized , however, that in our case H is not a linear operator unless the equation 
of motion is linear (h^'"' = for / > 0). In the following we will also include the 
possibility that the coefficients /i'-'-' depend on time. We will see below that crucial 
features of the non-relativistic limit also persist for the relativistic second-order 
equation (||). 

At some time to we specify an ensemble of initial conditions for Xm = Xm(to) by 
a probability distribution exp — S'o[x°]- At any later time t the generating functional 
for the n-point functions is given by 

Z[j, t]= f Dx' exp{-5o[/] + frnXn.it; /) + JmX*mit; x')} (8) 



with Xm{t', X°) the solution of the equation of motion (|I|) with the particular initial 
condition x(to) = X°- ^^ particular, the two-point function reads 

p™„(t) =< Xmit)xUt) >= Z~'—^\,=o (9) 

with Z = Z[j = 0] independent of t. More generally, for a given time t the state of 
the system can be characterized by the coefficients z^^'^^ (t) of a Taylor expansion of 

Z[j,t] 

OO OO -| 



with 

< x.,(t)...x,,(t)x;(t)...x;(t) >= ^-'4t'l.,....,(t) (11) 

The time dependence of Z[j,t] is determined by an exact evolution equation [?] 

dtZ[j] = frrrFrn[-^]Z[j] + J^F:^[^]Z[j] 

= -^ryh^^' AAA AAzM 

V" ^(0 Odd "^ "^ yr 1 

with dt the time derivative at fixed j. This is equivalent to an infinite system 
of evolution equations for the correlation functions z^^''), which obtains by taking 



appropriate partial derivatives of eq. (|T^ with respect to j. For the equation of 
motion (|I|) Fm involves at least one derivative with respect to j. In consequence, the 
evolution equation for z^'''''' involves the correlation functions z*-'^ '' -* with k'+l' > k+l, 
but not those with k'+l' < k+l. For a given k+l the correlation functions z^'' •' ^ with 
k' + l' > k + l act as a sort of "environment" for the evolution of z^'^''"': They influence 
the dynamics of z^'''''\ but are not affected themselves by the time evolution of z^'^''). 
Symmetries of Z[j] are conserved f\ by the evolution equation (|12]) if U is invariant 
under these symmetries. If Z involves initially only terms with an equal number of 
jm and j^, this property remains conserved by the evolution. More generally, the 
correlation functions z^'''^^ can be classified by their charge Q = k — I. Sectors with 
different charge are not mixed by the time evolution. If 5*0 [x°] in eq. (§) is neutral, 
only the neutral sector of Z[j] is nonvanishing. 

For conservative equations of motion (H), eq. (|T2D has a fixed point {dtZ^, = 
0) which corresponds to thermodynamic equilibrium at temperature T = jS"^ [Q, 
namely 

Z,[j,(3] = J Dxexpi-f3U[x] +fmXm+Jn.X*J (13) 

In this note we address the question if this fixed point is approached for t ^ oo. We 
will see that for generic initial conditions this is not the case due to the existence 
of infinitely many conserved quantities which characterize the evolution of the cor- 
relation functions. For the equation of motion (|T]) a first set of conserved quantities 
is easily found by observing that due to the factor i one has for arbitrary j 



- im§-H..F4])g£-zy\ (14) 

For J = one infers immediately that the trace of the two-point function (^ is 
conserved, dTrp/dt = 0. We also conclude that the operators (d'^/dj^djm) and 
iJpFp[-^] + jpF*[jr]) commute. This leads to an infinite set of conserved correlation 
functions Xjy 

q^^dz-^i-J^Yz] .0 (15) 



^More precisely, this holds only for symmetry transformations that do not mix j,„ with j* 
The transformations oi jm and Xm are related by the invariance oi j^Xm + jmXm- 



It follows immediately that the equilibrium fixed point ( [TB| ) cannot be reached un- 
less the initial values of all X^ exactly correspond to the ones of the equilibrium 
distribution! 

Actually, the X^ are only a small subset of the conserved correlation functions. 
We will concentrate in the following on the two-point function pmn and show that 
there are states for which this will behave exactly as the density matrix in quantum 
statistics. The evolution equation for the two-point function reads 



We note that only those terms in the general expansion of Z[j] (p!0|) which have zero 
charge 
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contribute to eq. (|16D 
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, , Pmn ^ -^ mn i iJ- rnn 
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Y - Z-^H \^^ ^'^ 

^ mn ^ ^-'■mpy 



dj dj;djn^.^^ 



oo 

•'mpqiri...qiri^pnriqi...riqi 



'^mpPpn + -^ 2^ i^mp. 
1=1 

We will consider states which obey 

Y -Y^ = [H, p] - i{K, p}, H = H\ K = K^ (19) 

This holds, in particular, whenever p is invertible, H — iK = Yp~^. In terms of the 
matrices H and K we can write 

^^p = ^^[H,p]-{K,p} (20) 

with Tr{Kp) = 0. The particular subclass of states which, in addition to the 
condition (|19|), obeys for all t 

{K,p} = Q (21) 

will be called quantum states . The time evolution of quantum states is governed by 
a von Neumann- type equation dp/dt = —i[H,p]. In particular, for linear equations 



of motion only h^'^^ contributes to Y ( ]TE| ) and all states are quantum states. For 
nonlinear equations of motion the Hamiltonian H still depends on the particular 
quantum state which is not only characterized by p, but also by higher correlation 
functions z^''\ Nevertheless, the commutator on the r.h.s. implies that the trace of 
arbitrary powers of p is conserved 

jTrip"^) = (22) 

This yields once more an infinite set of conserved quantities. Relation (^) is well 
known for the density matrix in quantum mechanics. For any function F{p) which 
admits a Taylor expansion it implies -^Tr F{p) = 0. In particular, the quantum 
mechanical entropy Sq = —Tr p In p is conserved. Linear quantum states are those 
for which H does not depend on p. 

The analogy of p with the quantum mechanical density matrix can be pushed 
even further. By the defintion (j^) p is hermitean and obeys the positivity conditions 
^mPmn^n — ^ ^^^ arbitrary complex vectors ^. In particular, the diagonal elements of 
p are positive semidefinite, p„„ > 0, and obey PmmPnn > |PmnP for all m and n. The 
expectation value of an arbitrary quantity which is bilinear in x, ^ = Xm^mnXn, 
can be written in terms of an operator associated to the matrix A^n as < ^ >= 
Tr{Ap). Different operators A,B do not necessarily commute. By suitable unitary 
transformations one can bring p into diagonal form p = diag{pn) with real positive 
semidefinite eigenvalues Pn > 0. Therefore p is invertible if none of the pn vanishes. 

In order to gain some intuition how restrictive is the "quantum condition" ( pi]) 
we consider an equation of motion which contains only quartic nonlinearities, i.e. 
/iW = for / > 2. For invertible p the conditon {K, p} = for a quantum state 
is equivalent to /C = or pY = Y'^p, as can be seen easily in a basis where p is 
diagonal. Then quantum states obey a condition for the four-point function z^'^^ 

7,(1) (o z(2) - n r(2) ^ - n (0^) 

Since this condition should hold for all t, the time derivative of the l.h.s. must 
also vanish. In turn, this gives a condition for z^^\ and so forth for higher time 
derivatives. Generically, quantum states therefore correspond to restrictions for all 
z^') . Nevertheless, z^gr has four indices and eq. (^) is an equation for a combination 
forming a matrix with two indices. One therefore expects that there are many 
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solutions of eq. (^) for zS'^'' with given h and p. The occurrence of quantum states 
obeying the conditions ([T9|) and ( PD seems quite generic ![] 

As an example for our discussion of quantum states, let us take the field equation 
(0) with a local quartic interaction W[^p) = |Ap^. In a basis where the index m 
denotes space coordinates (i.e. lattice sites), the nonlinear part of H takes the form 
^mligr — i^j'^mgl^mn^mr + SgnSgr)- The evolutiou cquatiou for the two-poiut function 
reads 

+^(< xix)x*ix)xix)x*iy) >-< xix)x*iy)xiy)x*iy) >)} (24) 

The conditions for a quantum state require that for every t there exists a function 
G{x,z) = G*{z,x) such that 

< x{x)x*{y)x{x)x*{x) > - < x{x)x*{y)x{y)x*{y) > 

dz{G{x, z) < x{z)x*{y) > -Giz, y) < xix)x*iz) >} (25) 

We will next consider a particular class of states where all correlation functions 
z^''^ are charge conjugation- and translation-invariant for / > 2. For these states the 
"quantum condition" ( p5D will be fulfilled by virtue of the symmetries with G{x, z) = 
0. Since U = ^ J d^x{diX*ix)dixix) + ^ix*ix)xix)Y} is translation-invariant, an 
initial translation invariance of the correlation functions i^'^ for / > 2 is preserved 
by the evolution. Due to the particular structure of the evolution equation (|T^), 
this property even holds if z*-^-* is not translation-invariant. (A translation- invariant 
"environment" for i*-^-* is not affected by the evolution of z^^' .) In a momentum basis 
the translation invariance of z^'^^ implies 

< X{Pi)X*{P2)x{P3)X*{Pa) >= A4{pi,P2,P3,PA){'2T^yS{pi-p2+P3-pA) (26) 

and the evolution equation for p{p,q) =< x{p)x*{l) > reads 

—p{p, q) = -^{P^ - (f)p{.P. l) (27) 

iX r dl^r d^s 



J (2'kY f27r')^ ^^^'^^' P' ^' ^ + ^ ~ P) - Ai{p,r,r + s - p, s)]{2TrY5{p - q) 



2M J {27rY{2'K 



''Mathematically, the simplest quantum state is z^-^^ ~ for I > 2. This seems, however, not a 
very natural state for physical systems. 



d f d'^s 



In addition, U is invariant under the charge conjugation x{p) ~^ X*{p)- Fo^' charge 
conjugation symmetric states one has ^4(^1,^2,^3,^4) = ^4(^2,^1,^4,^3) and the 
nonhnear term ~ A in eq. (^) vanishes. For a translation- and charge conjugation- 
invariant "environment" the evolution of the "density matrix" p is therefore governed 
by a linear von Neumann equation 
d , , . f d'^s 

W) 

with a hermitean free Hamiltonian H{j), q) = ■^{2t[Y5{jp — q). Despite the nonlin- 
earity of the microscopic equation of motion (|^) we have found a class of probability 
distributions for which the two-point function obeys a linear evolution equation! 
We should emphasize, however, that charge conjugation invariance of the initial dis- 
tribution does not guarantee that the state (or the environment for z^^^) remains 
invariant at later times. Charge conjugation maps x onto x* and is therefore not 
automatically conserved by the evolution equation. In general, only a subclass of 
states preserves charge conjugation symmetry. In particular, the thermal equilib- 



rium states (13) are translation- and charge conjugation-invariant. If we consider 
deviations from these equilibrium states where only the two-point function z^^^ dif- 
fers from the equilibrium value, the time evolution of p is governed by the von 
Neumann equation (|26|) . The superposition principle is valid, and for two solutions 
Pi and P2 the sum pi -|- p2 is also a solution. 

The construction of the general solution of the von Neumann equation (with 
arbitrary hermitean H) is well known from quantum mechanics: It is based on a set 
of wave functions -ip^^^t) which obey the Schrodinger equation 

^^i^)[t) = -iHi;^''\t) (29) 

with 

p{t)=Y.Pr.ij^-\t)®,lj^->{t) (30) 

Here p„ are the positive semidefinite time-independent eigenvalues of p and -ip^^^t) 
the associated eigenvectors in a standard normalization tp^^^*-ip^^^ = 5^"". (In a 
given basis for the degrees of freedom Xm we can write ■ip^'^\t) = anm{t)'4'm with 
time-independent basis vectors ipm such that Pmnif) = J^sPsO'smKn-) -^ pure linear 
quantum state obeys p^ = p Tr p. So far we have not chosen a particular normal- 
ization of p - the standard normalization Tr p = 1 can always be achieved by a 
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suitable rescaling of the degrees of freedom Xm- In summary, we have found a class of 
states for which the two-point function obeys all the dynamical rules of a quantum- 
mechanical density matrix. The time evolution of such states will not approach the 
thermodynamic equilibrium state. Just like in quantum mechanics, thermodynamic 
equilibrium can only be attained for suitable subsystems or macroscopic degrees of 
freedom. 

For the non-relativistic quartic scalar theory, we find from eq. (p8D that a trans- 
lation symmetric two-point function 

p{p,q) = A,{p){27rY6{p-q) (31) 

is conserved for arbitrary (positive semidefinite) A2{p) if the higher ra-point functions 
are the ones of a linear quantum state. We have therefore found an infinity of 
new fixed points of the evolution equation ([T^) . They are related to the thermal 
equihbrium state (|T3|) by 



Zi^'^^j] = Z^[j,f3] + zJ0^,r{p){A2ip) - Ai\pMp) (32) 

with A2J (p) the two-point function in thermal equilibrium. We emphasize that 
the connected higher n-point functions generated by W^ ^ [j] = In Zi ^ [j] differ 
from the ones in thermal equilibrium. (The disconnected parts are modified by the 
difference between A2 and A2J which implies that the connected parts must also 
be modified.) Most likely, the fixed points (^) are just a small subset of all the 
fixed points of eq. (0) - the top of an iceberg. Charge conjugation and translation 
symmetric states are not the most general quantum states. Furthermore, we have 
concentrated on two-point functions and not exploited the conserved quantities in 
higher n-point functions (cf. eq. (|1^)). This raises the question what singles out 
the thermodynamic fixed point - an issue which is probably related closely to the 
concept of entropy for suitable subsystems. 

Finally, one may wonder if the additional fixed points and conserved correlations 
appear only in the non-relativistic approximation (|^, or if they are also present for 
relativistic field equations. The second-order equation (^ can always be written as 
a coupled set of first order equations^ 

SSci 



ip{x) = 7r(x), 7r(x) 



6{p*{x) 



^In this form they arc not of the type (hi) with hermitean H. 



S,i = J d''x{diip*{x)dMx) + M'p{x) + -Xp\x)} (33) 

Consider first the particular two-point function^ 

pa{x, y) = ^ < '^ix)f*iy) - fix)'r^*iy) > (34) 

It obeys the evolution equation 

dtpA{x,y) = -i{-—{A^-Ay)pA{x,y) (35) 

+ ^(< V{x)f*{x)!f{x)(f*{y) > - < if{x)if*{y)if{y)if*{y) >)} 

and we immediately recognize the close similarity with eq. (0). Again, if the higher 
correlation functions are translation- and charge conjugation-symmetric, the term 
~ A vanishes and pA evolves according to the von Neumann equation. We empha- 
size that in contrast to the non-relativistic equation of motion charge conjugation 
symmetry remains conserved for the relativistic equation. A charge conjugation 
symmetric initial state always preserves this property at later time. (In this con- 
text it should be noted that the charge conjugation ip{p) -^ ip*{p),n{p) -^ tt*{p) 
discussed here is not equivalent in the non-relativistic limit to the transformation 
x{p) ~^ X*{p) discussed previously.) In conclusion, as far as pA is concerned all 
initial states with translation and charge conjugation symmetric 5*0 are linear quan- 
tum states. For these states Trp^ is conserved and the system cannot approach the 
equilibrium fixed point. 

The situation for the relativistic system is, however, more involved than for the 
non-relativistic approximation. This is seen from the time evolution of the other 
two-point functions 

p^ =< <^{x)<^*{y) >, P^ = J^ < ^i^)'^*{y) > 
Ps = ^ < n{x)^*{y) + ^{x)n*{y) > (36) 

which obey 

dtPv{x,y) = 2Mps{x,y) 

dtps{x,y) = M{p^{x,y)- p^{x,y)) + —{A^ + Ay)p^{x,y) 



^In terms of x(a;) (eq. (JSJ)) pA equals p up to corrections ^ x/M. 
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'< (p{x)ip*{x)ip{x)ip*{y) > + < ip{x)ip*{y)ip{y)ip*{y) >) 



2M 
dtpA^^y) = -'^Mps{x,y) + —{A^ + Ay)ps{x,y) 

+ ^i^x-^y)PA{x,y) (37) 

'< ip{x)ip*{x)ip{x)n*{y) > + < TT{x)ip*{y)ip{y)ip*{y) >) 



M2 



For A = the general solution is oscillatory wlieras for nonlinear systems the be- 
haviour of Pip,Pn and ps is not yet established. A partial approach to a fixed 
point is not excluded. In any case, we observe that for translation-symmetric 
Pa{p,<i) = ^2(p)(2vr)'^(5(p — q) the evolution equations for Pip^p-,, and ps are inde- 
pendent of Pa- Independence of pA holds trivially also for the higher correlation 
functions. Comparing with thermodynamic equilibrium where pA vanishes, we ob- 
tain again an infinite set of additional fixed points 

^i^^^j] = Z.[3.P]-^MZ j -^^A,{p){3:{p)3^{p)-3*M3.{P)) (38) 

where j*(p) and j*(p) are the sources conjugated to -n{p) and (/^(p), respectively. 

In conclusion, we find that there is no simple approach of conservative classical 
systems to thermodynamic equilibrium. In many respects the analogy of the evolu- 
tion of classical correlation functions with quantum mechanics is striking. We are 
all used to the picture that classical mechanics arises from quantum mechanics in 
certain limiting cases. Our findings suggest that the opposite may also be true: The 
dynamics of quantum mechanics and perhaps also quantum statistics follow from a 
statistical treatment of classical systems in particular ensembles! 
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